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I. INTRODUCTION 


Studies in the western North Pacific (Jarrell) ct Ն M 
1978) and the eastern North Pacific (Thompson, et al... oe 
link short-term tropical cyclone track forecast errors, in 
part, to an incorrect initial position. Prior to requesting 
any of the various objective forecast aids (persistence, 
climatology, analog, statistic or dynamic), a current 
"warning position” must be determined. Initial positioning 
error is defined as the distance between an initial warning 
position and the corresponding post-season Best Track (BT) 
position. Neumann and Pelissier (1981) eliminated the 
effect of this initial position error in their evaluation of 
forecast errors by shifting the forecast tracks of the 
various objective aids so that the initial position was 
coincident with the BT position. 

Efforts continue to be made to reduce this initial posi- 
tionineg terror Prior to 1983, the Joint Typhoon Warning 
Center (JTWC), Guam, used the speed and movement based on 
the latest fix information to extrapolate a 6-h warning 
position. This extrapolated position was used as the 
"current warning position” when requesting forecast aids. 
The extrapolation procedure introduces uncertainty in this 
critical short-term forecast variable. Since 1983, the JTWC 
minimizes this uncertainty by using the most recent "working 
best track" position as the initial warning position for the 
purpose of requesting objective aids. The "working best 
track" (WBT) is the adjustment of the previous warning posi- 


tions during each warning cycle to better approximate the 


recent storm track based on the accumulated information. In 
most cases, positions 36 to 48 hours prior to the most 
recent position are not ‘Significantly altered. This 
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procedure is similar to the approach in determining the 
Semtieial best track’ after completion of the storm. WBT 
positions at -6 h, -12 h, -18 h, and -30 h are sent to the 
Fleet Numerical Oceanography Center (FNOC) to generate the 
various climatological and dynamic forecast aids that are 
used in preparing the warning position and future warning 
track. The +6-h forecast positions from the objective aids 
are blended with the most recent fix positions to determine 
mie. Current warning position’. Greater confidence should 
result since extrapolated tracks from each fix are not the 
only information utilized to establish the warning position. 
Rather, the forecast tracks from the various aids provide a 
kind of hindsight for judging the likely accuracy of each 
Ex, Subjectivity still enters the procedure during the 
blending process. The Typhoon Duty Officer (TDO) must 
consider the likely accuracy of the fixes from the various 
platforms (satellite, aircraft, radar or synoptic) in deter- 
mining the optimal warning position based on his experience 
and recent platform performance, i.e. which’ platform has 
given the best indication of recent storm movement. The TDO 
has available one to 15 fix positions during the six hours 
Since the previous warning. 

The goal of this thesis is to develop an objective 
procedure for the TDO to use in determining the initial 
PFacnıng position. The strategy for the objective warning 
position determination proposed here is based on the idea 
that it is often easier to determine which of the storm 


center fixes to accept if the forecaster knows the future 


Erack. That is, hindsight often allows the forecaster to 
select more intelligently from a number of possible fix 
positions. The TDO might then use the warning position from. 
the objective scheme as a "first-guess" position. This 


position is then adjusted to reflect consistency with the 


synoptic reasoning that forms the basis for the forecast 
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track. It is hoped that this objective warning position 
could be used as a "tool" by the tropical cyclone ‘forecaster 
to provide a more consistent initial warning position and a 
more accurate short-term forecast. 

In this thesis, a detailed description of the objective 
technique procedures and weighting functions applied to the 
individual fixes will first be given, followed by a summary 
of the sensitivity tests applied to the dependent data 
sample. Results from both dependent and independent samples 


will then be discussed. 


p 


I ROCEDURES TOFT THE OBJECTIVE TECHNIQUE 


Objective initial positioning has been proposed by 
Morford (1979) as an essential step in improving the short- 
term (less than 24-h) forecasts at JTWC. Simpson (1971) had 
earlier proposed a decision-tree format for establishing the 
EWtral warning position. This decision-tree approach 
attempted to use objective aids and similar reasoning in 
each forecast cycle to insure consistency. However, 
Simpson's technique was subjective rather than objective. 
The objective scheme proposed here simulates hindsight by 
estimating the future positions associated with each fix 
through an economical and viable short-term forecast tech- 
nique, the western North Pacific CLIPER, described below. 

The western North Pacific GEFPER (CLimatology .and 
PERsistence), which was developed by Xue and Neumann (1984) 
at the National Hurricane Center (NHC), uses regression 
equations to relate future storm displacement (DISP) to 


eight basic environmental predictors: 


DISP = f( Xo, Yo, pa 2 27 NM ; n" ՐԹ") 
Xo - Initial longitude 

Yo - Initial latitude 

,»։».-. 1.7 lO-west translation 
۳۶۰۰٠٦٦۰٠0 ٦73 beast too west Exanslation 
Eee ከ co north translation 
2۰۰۰۰۶ 3۰٠٠٠٠٠٦ Mh SOUL tos north translation 
W = Initial maximum wind speed (kt) 

D = Julian date. 


Note that westward and northward translations are defined to 
be positive as the typical cyclone track is toward the 


northwest in the western North Pacific. 
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Because storm track forecasts to 24 h rely heavily օո 
persistence, the primary input to CLIPER is the past 12 and 
24 hour movement. However, higher order (up to third- 
order) terms can also serve as predictors in the regression 
equations. when all eight basic parameters are included, 
there are 165 possible products and cross-products as poten- 
tial predictors in the regression equations. The applicable 
regression coefficients were derived for storms south of 
35°N and west of 150°E during the months of May through 


December (see Xue and Neumann, 1984, for details). 


A.  INTERPOLATION TECHNIQUE 


Upon receipt of a fix position, past positions at fix- 
time minus 12 hours (F-12) and fix-time minus 24 (F-24) 
hours must be derived to provide the required input to 
CLIPER (See Figs. 2? P: Rather than linearly interpolating 
between the -12 and -24 h WBT positions, a smooth WBT is 
determined from which the desired positions (F-12, F-24) can 
be interpolated. A third-order polynomial is fitted to the 
latest warning position and the -6, -12, -18 and -24 h WBT 
positions. As stated above, the TDO will adjust the WBT as 
far back as 36 to 48 hours to better reflect the latest 
information. As the greatest confidence can be placed on 
the "earlier" positions, higher weighting factors are given 
to the -12, -18 and -24 h positions (Table 1). The polyno- 
mial coefficients derived from this fitting routine allow 
the determination of interpolated positions at any time 
along the smooth WBT. For example, if fix A (Fig. 2.1) was 
received at 0430 GMT, position B at 1630 GMT (F-12) and C at 
0430 GMT (F-24) would be interpolated using the third-order 
polynomial coefficients. These three positions, the current 
maximum wind speed (taken from the most recent warning or 


from the fix if it is from aircraft reconnaissance) and the 
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mre Julian date/time are input to CLIPER to generate a 
72-h forecast. Since the. goal is an improved initial 
warning position, only the 12-h and 24-h forecast positions 
are calculated. However, the 36-h through /2-h track is 
available if the complete forecast is desired. Linear 
interpolation is used to derive CLIPER forecast positions at 
۰51110 +18 h from fix time. 


TABLE 1 
5... 8 070 
Time (Hours) 
-24 -18 -12 -06 00 
Weight: 20 15 ՆՏ 10 10 


ባኮ]! order polynomial fit of the future  CLIPER 
warning positions (W+6 to W+24) and prior (00 to W-24) WBT 
positions is used to determine a smooth estimate of the 
storm movement. The polynomial routine allows for user- 
specified weights at each fitted position. Different order 
polynomial fitting will be included as an option to be 
selected by the TDO in the interactive version of the 
scheme. Larger weighting factors (Table 2) are given to the 
prior positions to assure a smooth evolution from these 
relatively well-known positions. Separate polynomial 
fitting of the latitude and longitude positions with time 
was adopted as an alternative to fitting the time sequence 
of latitude/longitude pairs. These time-dependent, fourth- 
order polynomial coefficients are used to determine the 
tentative warning position (the +6-h position in Fig. 2.1). 
Notice that the fix position (A) is not included in the 
polynomial routine. If the fix is close to the previous 


warning time or an erroneous fix is encountered, a large 


5ا 
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forecast and interpolated positions. 


movement discontinuity might be inferred and could not be 


fitted by a fourth-order polynomial. 


TABLE 2 
Weighting factors applied when determining a tentative 
warning position from each fix. I is the time 


increment between the last warning and the fix. 
Time (Hours) 


-24 --189 0232 06 +941 ٦٦ 
Weight: 20 50 40 10 10 5 1 1 i 
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when more than one fix is available, a weighted average 
of the interpolated positions (position D in Fig. 2.2) gives 
the first iteration of the warning position. The procedure 
for determining the TENES given to different types of 
fixes will be described in Section 2Շ. Up to 10 fixes may 


be included for each warning position determination. 


Pa dieo E 135 E 140° E 145” E 





RAS 2.2 Determination of arang M tren (Dy with 
ነ! ”  ረ፡ዝ፡ቋ፡ርፎ available. BERG „are 1 > 
positions; 'o', working best track. 


Consideration is also given to the potential impact of 
positioning errors of the fix platforms. In a second itera- 
tion, four adjacent positions are generated from the first 
Meera eion warning position (point D in Fig. 2.2) by adding 


1 


an "observational error’ in each Cardinal direction (Figs. 
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2.3. and Անել. Each error position is treated as a "fix" 
position with CLIPER forecasts generated and a polynomial 
curve fit as described above. The weights applied at the 
specified times in the polynomial curve for the second iter- 
ation warning position were changed due to the increase in 
the number of fitted points (Table 3). The four new tenta- 
tive positions are weighted equally since each would have an 
equal probability. The final warning position is determined 
as the arithmetic average of these four new positions (Fig. 
2.3b). Sensitivity tests show that smaller objective warning 
position errors from the best track result if the observa- 
tional error iteration is incorporated. Essentially the 
same improvement occurs if the observational error is 10,15 
or 20n mi. Therefore, an observational error of 15 n mi is 


taken as the default value. 


14° N 


125°E 130 E 135” E 140” E 145 E 





ace 2.3a Observational error positions (labeled 1-4) 
added in cardinal directions and the previous storm 
positions B and C at -12 and -24 hours. 
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Uc 2535 Second iteration of the objective warning 
Posttion 'e 


TABLE 3 
Weighting factors when fitting the second 
lteration warning position 
Time (Hours) 
-24 -18 -12 -06 00 +06 +12 +18 +24 +30 
Weight: 25 75 45 45 20 իջ 5 1 1 1 


B. INITIAL WARNING PROCEDURE 


A special procedure is necessary to start the objective 
technique with a new storm since the prior "working best 
track" positions required as input to CLIPER are not avail- 
able. To generate CLIPER forecasts for the first four warn- 
ings, a 24-h history must be developed. All that is 
required in the objective procedure is a fix about 6 h prior 


to the first warning time, e.g. a synoptic fix determined by 
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the TDO. A rhumbline course and speed between the most 
recent fix and the -6 h estimated position is used to extra- 
polate backward 24 h for the required history. As the 
warning sequence continues, WBT positions can be included as 


the -6, -12, -18, or -24 h positions become available. 


C. WEIGHTING FACTORS 


In the process of blending the various fix positions, 
the TDO subjectively weights each fix based on the fix 
platform accuracy (i.e. aircraft fixes are generally more 
accurate than satellite, e.g. Jarrell,1978) and the fix 
time. More recent fixes are normally given greater 
emphasis. 

Weighting factors for the aircraft, satellite, radar or 
synoptic fixes àre required for the objective technique to 
represent the expected accuracy of each fix (Table 4). The 
satellite analyst assigns an accuracy estimate (Position 
Confidence Number, PCN) based on the tropical cyclone cloud 
signature. A PCN value of 1 would be the most accurate 
estimate while PCN 6 would have the least accuracy. 
Aircraft fixes are grouped according to a combination al 
navigation accuracy and meteorological accuracy. The latter 
is a largely subjective estimate by the Air Reconnaissance 


Weather Officer (ARWO) of the accuracy of the wind/pressure 


center location. No indication of the accuracy of the radar 
fixes is provided, so that all radar fixes are grouped in a 
single category. Due to the very small number of synoptic 


fixes, they are grouped with the radar fixes for the purpose 
of determining accuracy. 

JTWC annually evaluates the mean accuracy of satellite, 
aircraft reconnaissance and radar fixes by comparing them 
with the BT position at the corresponding time (JTWC Annual 
Typhoon Reports, 1981-1983). A Similar approach has been 
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TABLE 4 


Grouping of the fix platforms for determining 
accuracy from 19812 - 1953 fixes 


Satellite: PCN 5 and 6 loose organization) 
PCN 3 and 4 well defined organization) 
PEN 12514222 eye present) 


Aircraft: Group Navigation Meteorological 
accuracy . accuracy. 
<7 <9 <7 nmi 268 ۱9 7 
<7>9 < 7 > 
>7<9 = 7 < 9 
»7»9 » 7 » 9 


Radar and Synoptic 


used with all storms during 1981- 1983 to extract the 
medians and standard deviations as potential weighting 
factors for each fix type. To determine the fix accuracy, 
the time difference between the last known position and the 
latest fix (hours and minutes) is converted to a percentage 
of the 6-h increment. This percentage is then used to 
linearly interpolate the point along the BT that corresponds 
to the fix time. The distance between the fix and this 
corresponding BT position then determines the fix accuracy. 
Table 5 is a list of the means, medians and standard 
deviations for the 1981-1983 storm seasons. Ill-defined 
storms, multiple centers,  upper-lower layer cloud signature 
decoupling, etc., in a small percentage of fixes contribute 
to large displacements from the official BT. These outliers 


serve to shift the distribution and bias both the means and 


the standard deviations. Thus, the median in each group is 
chosen as a more satisfactory measure of accuracy. As shown 
in Table 5, the distinction between the 3-yr average median 


ከ ከ he most accurate fix platform (Aircraft <7 <9 ) and 
least accurate (Satellite, PCN 5 & 6) is only a factor of 
three. To provide more discrimination between fix types, a 


third power of the median is utilized. 
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Consider an example of how the fix platform weighting 
factors might be applied. Prior to warning time, two satel- 
lite fixes (PCN 1 & 6 ) and one aircraft fix "7 ane 
received. The determination of the weighted mean position 
(Fig. 2.2) is then given by the sum over i fixes of: 

X= > (XI "Mi ` `. 


and 
ፕ = 211925 MI A 
where Xi and Yi are the latitude and longitude associated 
with each fix. Mi is the reciprocal of the median accuracy 
for the fix type i raised to the third power. In the 
example: PCN = 1, Mi = 1/8.5 and Mi = 4y614.1; V PONE 
Mi = 1/18.8 and Mi’ = 1/6644.7;  Nav/Met accuracy = <7 29m 


Mi = 1/6.8 and Mi ՀԵՆ արիր ՏՍ 2 2E 
Warning position will be heavily slanted toward the aircraft 
and PCN = 1 satellite fixes. 

The second weighting factor takes into account the fix 
receipt time. The TDO normally places greater confidence on 
fixes obtained close to the desired warning time, especially 
if it agrees with previous expectations. Warning No. 29 for 
Typhoon Marge during 1983 (Fig. 2.4) illustrates the need 
for this time bias. Two satellite fixes of equal position 
confidence (PCN = 2) were received at 12 GMT and 16 GMT in 
support of the 18 GMT warning. Based on the official fares 
cast, the JTWC must have placed most of their emphasis on 
the later fix. The objective routine without a time 
weighting factor would equally weigh the two satellite 
fixes, which would result in a much slower movement. 
Therefore, a linear weighting function is assumed using the 
time difference between the most recent warning position and 
the fix time (Fig. 2.5). Fix positions obtained five to six 
hours after the last warning time will be given ten times 
more weight than fix information received within the first 


hour after that warning. A second-order weighting function 
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was tested, but the improvement was not significant. Two 
weighting factors (fix platform and time of receipt) are 
then included in the determination of the weighted mean 
position for i fixes as follows: 
X > (Xi = mM) Ի ԱՋՆ Ժ 0 ՈւՄ 
and 
Y > (Yi * Mi)  + +187۳7 


where Mi is again the reciprocal of the median displacement 


raised to the third power and Ti is the time bias. 


D. SENSITIVITY TESTS 


The goal of the present technique is to synthesize 
objectively the TDO's procedure in deriving a consistent and 
accurate warning position. The skill of the objective tech- 
nique is determined by comparing the objective warning posi- 
tion and JTWC warning position to the corresponding BT 
position. The ultimate goal is to provide a warning posi- 
tion that will result in reduced Short-term forecast ۶۶'٣٦ 
from the objective aids (in our case, 000 ٦ 

Typhoon Pat (1982) was selected as the initial test case 
because it had large warning and forecast errors (Fig. 2.6). 
The storm tracked to the west-southwest while still in the 
tropics, decelerated during the recurvature stage and then 
accelerated following recurvature. It can be seen that 
Typhoon Pat's movement departs significantly from a smooth, 
uniform track. Seven additional storms (Thad and Bill- 
1981; Dot and Gordon - 1982; Marge, Herbert and Abby - 1983) 
were added to the test data base for a total of 226 warning 


positions. 


Լ. Extrapolation versus Interpolation 


An extrapolation approach to applying the CLIPER 


technique to each fix was also tested. A tentative 6-h 
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٦۷٦۷۹ 9 E phoon Marge (1982) 


een un errorem © ective „varning Position 

Es Go exe ltsi10on of a, time las. , bes rack 

xs fir position; '*', objective warning position; 
x JTWC warning position 


warning Es is determined based on a rhumbline course 
and speed (dashed line in Fig. 2.7) between the fix and the 
ծ... թո 12-h prior to the desired warning position. 
Using the prior 12-h rather than the 6-h position minimizes 


eea extrapolation angles that might occur if the fix 
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WEIGHT 


TIME (hr) 





Fier 0 Weighting factors for fixes as a function 
of the time since the previous warning. 


receipt time is early in the 6-h interval. As the required 


CLIPER inputs are the current, -12 and -24 h positions, the 


corresponding values in Fig. 2.7 would be the extrapolated 
warning positions (A or B) and the 18 GMT and the 06 GMT 
positions from the best track. No time interpolation along 


the best track is necessary in this approach, in contrast to 


the interpolation method shown in Fig. 2.1. Although this 
extrapolation approach is easy to apply, it exaggerates 
erroneous fixes by the extrapolation. The interpolation 


method described in Section 2.2 is more conservative as the 
desired -12 and -24 h positions are interpolated from rela- 


tively well known positions along the WBT. 


2. Polynomial Curve Fitting 


Both third- and fourth-order polynomial curve 
fitting routines were evaluated in determining the first and 
second iteration objective warning  positons (see Figs. 2.2 
and 2.3). Experimentation with the eight-storm data base 
indicated that choosing the fourth-order polynomial resulted 
in the smallest initial warning errors. Third-order polyno- 
mial curves did not approximate the working BT and forecast 
track as well and resulted in larger average warning 


position errors in most cases (Table ር Higher order 
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polynomials would require additional positions to 
the smooth curve and tend to add complication 


significant improvement. 
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06002 based on fixes at 04002 (A) and at 01302 (B). 


3. Weighting Functions 


The next step in the sensitivity analysis concerned 
a “Goodness or Fir” approach. If the smooth polynomial 
curve correctly fits the past and forecast positions, it is 


presumed that a better approximation to the official BT will 


follow: Greater emphasis is given to the prior (-6 to -24 
h) positions since they are considered "known" at warning 
time. The evaluation was to determine whether larger 


weights should be given to the later forecast positions (+18 


and +24 hour), to the earliest forecast positions (warning 
position and +6 hour), or to treat equally all forecast 
positions. Larger weights at the tentative warning position 
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TABLE 6 


Average warning position error (n mi) 
with 3rd or 4th order polynomial curves 


Storm No. of JTWC Polynomial Order 
Warnings Jd 4th 
TY Pat en 24 24.9 1973 18.4 
STY Abby § 3) 50 131٦ 16.0 15.0 
TY Thad ( 28 296 2608 26.1 
TY Marge (8 4 30 18.4 24.1 23.1 
TY Gordon (8 | 37 15.5 16.3 15.4 
TS Herbert (83) 10 14.3 12.6 127 
TY Bill (81 16 18.7 Ike 2 17.3 
EE Dot (82) 31 | 2099 18.7 
Total 226 
Weighted Average 17.6 19.4 18.4 


and the +6-h forecast generate better results in six out of 
the eight storms in the test data base (Table 7). This 
increase in accuracy of initial warning position is consis- 
tent with the expected importance of persistence. 

Table 8 is a list of the average warning position 
errors after inclusion of the time bias (Section 2C). The 
early forecast position weighting emphasis in Table 7 was 
adopted for these tests. All eight storms had reduced 
warning position errors when the time weighting factor in 


Fig. 2.5 was included. 


29 


TABLE 


Average warning position error (n.mi. 
data base emphasizing early (00/+06) 


Of che 0 ص9ج‎ 
+ 
warning posıtions 


) 
or 1ጋቢሮክስሽስ 724) 


Storm No. of 


JTWC Weighting Emphasis 

Warnings ar 06 

TY Pat on 4 24. 18.6 1075 

Y Abby | 3) 50 Lig 1291 13.6 

Thad 11 28 2215 265 278 

TY Marge (8 ) 30 18.4 239 22.4 

TY Gordon (8 } 37 15.5 116 14.8 

TS Herbert (83) 10 14.3 12085 14.8 

TY Bill ի 16 1887 17.4 15.9 

TY Dot (82) d 170 190 18.5 
Total 220 

Weighted Average 176 18 7 7 

TABLE 8 


Average warning AO errors (n mi) after 
inclusion of the time weighting factors 


Storm No. of JTWC Inclusion of 
Warnings Time Bias 
TY Pat on 4 24.9 18.0 
STY Abby | 3) 50 ried E7 
۳ئ‎ 2. HE 12 
arge 
IY Gordon (82) 37 3.5 14.0 
ր ՀԱՐՈ) ከ 18:7 15:5 
i : 
TY Dot (82) 1 17.1 16.8 
Total 226 


ԷԶ 
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Weighted Average 
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III. RESULTS AND DISCUSSION 


In this section, the objective technique will be 
compared to JTWC's procedure with respect to the average 
warning position error and the resulting 24-h forecast posi- 
tion. Several storms will be highlighted to illustrate the 
strengths and weaknesses of the objective scheme. A strati- 
fication of the storms by intensity will also be illus- 
trated. 

۱ Subsequent to the testing phase (described in Section 
Poe 22 additional storms were run for a total of 637 inde- 
pendent warning positions. Inclusion or exclusion of a 
particular storm from the three-year data base is based on 
the following constraints: 

(1) A fix position 6 h prior to the first warning 
(to initiate the objective technique) must 
exist; 

(2) Since the current objective technique is 
limited to 10 fixes per warning, storms with 
a large number of radar fixes per warning 
were not included; 

(3) Storms with a majority of their path outside 
the latitude/longitude domain of the CLIPER 
regression equations were not used; and 

(4) Several storms that were either short-lived 
or included periods of time during which 
JTWC warnings were unavailable were not 


included. 
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A. INITIAL POSITION ERRORS 


The accuracy and consistency of the objective technique 
as compared to JTWC can be measured by the means and stan- 
dard deviations of the initial position errors. Table 9 
lists these statistics for the 8-storm test base and indi- 
cates if the objective position error for each storm is 
smaller (W), the same as (T) or larger (L) than ٦ 
Tables 11 - 13 indicate the same statistics for the 22-storm 
independent sample. As shown in these summaries, the objec- 
tive technique produced a more accurate warning position in 
1981 and 1982. The objective technique resulted in smaller 
errors for 6 of 7 storms during 1981 and 3 ot 7 storms 
80176 ISA The percentage of storms in the win and tie 
categories represents 100% of the 1981 storms and 71% of the 
19982 “stormen Within the independent data set, Typhoon 
Nelson (1982) represents the best performance of the objec- 
tive technique while Typhoon Bess (1982) had the worst 
record compared to JTWC. 

Presumably, the change in the JIWC warning procedure 
during 1983 resulted in smaller overall errors compared to 
the objective scheme for the 1983 storms. This change may 
also account for the improvement in JTWC's performance in 
1983 compared to 1982 and 1981. À reduction in the number 
of storms in which the objective scheme produced smaller 
errors (3 of 8 storms) results in a win and tie percentage 
of only 501 Tor LOST Continued testing of the objective 
technique with 1984 storms is required to determine whether 
1983 was an anomalous year. 

The yearly summaries also indicate that the objective 
technique will produce a more consistent warning position, 
Since the standard deviations were smaller overall for all 
three years as compared to JTWC. Tables 11 - 13 indicate 
that 5 of 7 storms in 1981, 4 of 7 storms in 1982 and 4 of 8 
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[. SS 


Summary ot wasuinegsposition errors (ուտլ.) for 
JTWC JT) and gosse ive (5350; fOr Tod 8-storm 
test Se base. win առ 1) orsloss E 


for the Sure. d is M dea and 
the Student-t score is given. 


Storm Nom cot BU Std. Dev WTL t 
Warnings JT BW JIT OBJ 
ur Dot (82) 31 20 20 L3 14 1 ՀԱԱ 
TS Bert 9133) 10 18 ING 10 -5 W -. 44 
EY 48 50 11 13 10 9 Ն 81. 
ՍՄ mae (8 24 27 18 16 10 W -2.47* 
w Thad (8 à 28 26 27 J1 36 L 20 
TY Marge (8 4 30 20 22 17 18 E 27 
TY Gordon (82) 37 15 14 10 9 W -.67 
ዜ Bill (81) 16 19 15 13 12 W -.84 
Total 258 
Weighted Average 185 17.9, pa Bob ፌራ 


* Difference between the objective error and JTWC 
mac errom is sienificant 35059955 confidence 
eve 


storms in 1983 had smaller standard deviations. The 
Student-t test is made for each storm to test whether the 
difference between the objective warning position error and 
that of the JTWC is significant at the 95% confidence level. 
For the 22-storm independent data base, only Typhoon Nelson 
41982 ) would have had significant reductions in warning 
position error if the objective technique had been used. 
However, the warning positions of Tropical Storm Ben (1983), 
Super Typhoon Forrest (1983) and Tropical Storm Georgia 
(1983) would have been significantly degraded. Examples 
from these four storms and others of the 22-storm sample 
will now be examined to determine when the objective tech- 
nique should or should not be expected to provide accurate 
Bnortral position guidance. 

Deere 2.0 indicates the relationship of the official BT, 
the JTWC warning positions and the objective technique 
warning "positons for Typhoon Pat (1982). The objective 
technique performs very well during the recurvature (change 


of a dominant northwest to northeast movement around the 
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subtropical ridge axis) and seems to provide a very consis- 
tent acceleration track after recurvature. Ten of the 0 
storms analysed were recurving SyStens. Of these ten, six 
storms had overall smaller initial warning position errors 
when employing the objective technique. However, a problem 
does arise in situations where rapid acceleration follows a 
slow movement during recurvature.: For example, the. CLIPER 


forecast for the 18 GMT 22 August 1981 warning position of 


Typhoon Thad (Fig. 3.1) is based on the slow recurvature 
around the ridge axis from 12 GMT 21 August through 06 GMT 
22 August. During the post-season analysis, the TDO seems 


to have discounted the fix information (positions A,B and C 
in Fig. 3.1) to arrive at the accelerated 18 GMT 22 August 
warning position. It does appear that the objective warning 
technique based on CLIPER may be too conservative during 
these rapid acceleration cases. 


Typhoon Ben (1982) illustrates a scenario in which the 


CLIPER-based objective technique will not perform well. As 
can be seen in Fig. 3.2, the latitudes of the objective 
warning positions are quite accurate. However, the CLIPER 


model does not anticipate the non-climatological westward 
movement in the northern latitudes. 


Typhoon Nelson (1982) is characteristic of systems in 


which the objective technique Should produce accurate 
Warning positions. The first two-thirds of Typhoon Nelson's 
track, (ETS SN) is consistent with a persistence-type 
track. Objective aids which stress persistence for gaa. 


short-term forecast, such as CLIPER, are expected to verify 
well in these situations. Throughout this portion of 
Typhoon Nelson's track, the average JTWC warning position 
error was 21.4 n mi with a standard deviation of 14.1 n mi, 
while the objective scheme resulted in an average error of 
15.7 nmi with a standard deviation of 9.7 nmi The later 


portion of Nelson's 6۲407۰٦ BT included a 
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counterclockwise loop. An objective scheme based highly on 
persistence would not be expected to handle looping tracks 


very well. Due to the slow forward speed between warnings, 
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the objective technique still had a smaller average warning 
Position error (20.2 nni, seandard.deviatıon 11.5 n mi) 
Bonpared toWTWE" (average position error 26:1 n mi, standard 
feviation 20.2 n mi). Unfortunately, this skill in looping 
Situations does not hold for all of the cases. Six other 
storms (Pamela during 1982 and Percy, Lex, Abby, Bess and 
Sperry. during 1983) had. a loop sometime during their 
lifespan. Only three of the looping storms (Nelson, Sperry 
and Bess) had smaller warning position errors during the 
looping phase when using the objective scheme. In each of 
these three storms, there was relatively slow movement 
through the loop. In the remaining cases, the objective 
technique based on CLIPER did not handle the rapid direction 
changes associated with a tight loop. 

The use of the objective technique should not be ruled 
pur in all non-climatological situations. As an example, 
the entire lifespan of Tropical Storm Sperry (1983) 
consisted of a clockwise loop in the region east of the 
Philippines (Fig. Hed Once again the slow speeds 
throughout the majority of the loop allowed the objective 
technique to produce smaller warning position errors than 
JTWC (see Table 13). 

The TDO may have to adjust the objective warning posi- 
tion when additional knowledge not reflected in the fix 
positions is available. Such a situation appeared for the 
warning position of Typhoon Thad on 00 GMT 23 August 1981. 
The objective technique relied on two satellite fixes shown 
in Fig. 3.6 Although the resulting objective warning posi- 
tion was consistent with the prior storm track, the position 
error was 1057 n mi compared to a JTWC error of 30 nmi 
Evidently, the TDO had more information than the two satel- 
lite fixes provided. 

In summary, the 30-storm sample provides a fair repre- 


sentation of the various forecasting scenarios including 
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recurvature and looping cases as well as straight tracks. 
It has been shown that the objective technique will provide 
a reliable warning position in most situations, including 
some cases in which slow loops or erratic movement is 
experienced. In an operational version of the object 
technique, the TDO should have the capability of interac- 
tively  bogusing additional information for situations in 
which the fixes do not reflect the best estimate of the 


actual storm location. 


B. 24-H FORECAST EVALUATION 


A second comparison between the objective and JTWC 
warning position is to determine which method provides the 
most accurate 24-h forecasts. The JTWC and the objective 
warning positions were used with the -12 and -24 h best 
track positions in the CLIPER scheme to forecast the +24 h 
position. It was expected that initiating CLIPER with more 
accurate warning positions would result in more accurate 
short-term forecasts when compared to the +24 h BT position. 
The control case is a CLIPER forecast generated entirely 
with BT positions which should produce better short-term 
forecasts than either the JTWC or objective warning posi-: 
tions. Because the first 24 h of each storm track is 
required as past history in CLIPER, the first four warning 
positions in each sample cannot be evaluated. Furthermore, 
the warning positions during the last 24 h of each storm are 
eliminated to allow for the +24 h verification positions. 
The total number of verifiable warning positions in the 
sample was reduced to 456 and the short-lived Tropical Storm 
Sperry was eliminated in this evaluation. 

Tables 14 - 16 are listings of the 1981 - 1983 summaries 
in the same format as the tables for the initial position 


error evaluation. As expected, the 24-h forecasts based 
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entirely on BT positions were the most accurate. As in the 
initial position evaluation above, the 24-h forecasts based 
on the objective warning positions compared favorably with 
those based on JTWC warning positions during 1981 and 1982. 
The average 24-h forecast error for the objective scheme was 
۱۶۰۰ mi for 1981 and 113 ኳ mi for 1982. This can be 
compared to 127 n mi for 1981 and 111 n mi for 1982 based on 
ENNUG s" "warning "positions. The percentage of win and tie 
category storms was 55% in 1981 and 60% in 1982. However, 
only 20% of the storms during 1983 had better forecasts from 
the objective warning positions. According to the Student-t 
scores in Tables 14 - 16, none of the differences between 
the objective technique forecast errors and those of the 
JTWC were significant at the 95% confidence level. The 
standard deviations in Tables 14 - 16 indicate the 24-h 
forecasts based on the objective warning positions were 
slightly more erratic for 1982 (65 nmi) and 1983 (66 n mi) 
compared to forecasts from JTWC warning positions during 
1982 (62.9 n mi) and 1983 (60 n mi) 

The initial position error categories (Tables 11 - 13) 
are compared in Table 10 to the 24-h forecast position 
categories from Tables 14 - 16. Optimally, all storms in 
the win category of the initial position evaluation would be 
expected to result in a win category in the 24-h forecast 
position evaluation. It was not expected that storms with 
initial position losses would result in 24-h forecast posi- 
tion wins. As indicated in Table 10, most of the storms do 
fall in the win-win and loss-loss categories. However, two 
cases (Super Typhoon Marge (1983) and Typhoon Ken (1982)) 
fall in the loss-win category, i.e. an initial position 
error loss for the objective technique results in a 24-h 
forecast position win. Both storms were characterized by 
recurving tracks. A relatively large number of storms (6) 


with a win in terms of smaller initial position errors 
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nevertheless had a corresponding loss in terms of a 24-h 
forecast comparison. Of the six storms, three (Elsie 
(1981), Bill (1981) and Thelma (1983)) were recurving 
systems. Two others (Kit (1981) and Nelson (1982 )) were 


basically westward tracks, although Kit contained two 
radical course changes. Tip (1983) had a short-lived track 
through the South China Sea. Fig. 3.7 illustrates one 
possible explanation of the disparity. While the objective 


scheme produces a more accurate warning position (A in Fig. 
3.7) compared to JTWC, the recurvature to the northeast is 
forecast too late if the objective warning position is used 
in CLIPER. 

The discussion above does not support the expected 
coupling between the more accurate initial warnıng position 
and a more accurate short-term forecast position. This is 
illustrated in the following storms. The objective tech- 
nique resulted in significantly smaller initial warning 
position error (objective technique, 17 n mi; JTWC, 23 n mi) 
for Typhoon Nelson (see Table 12). However, the average 
24-h forecast error from the objective warning positions(124 
n mi) was larger than forecast errors based on the JTWC 
warning positions (111 n mi). Table 17 is a list of the 
CLIPER 24-h forecast errors for Typhoon Nelson based on the 
BT position, JTWC warning position and objective warning 
position. In this case, CLIPER forecasts with objective 
initial positions provide less accurate forecasts from 12 
GMT 24 March 1982 through 12 GMT 27 March 1982 (warnings 25 
- 36, Table 17) as Typhoon Nelson passes over the 
Philippines. The CLIPER forecasts from objective positions 
are again less accurate than from JTWC positions from 18 GMT 
29 March through 06 GMT 31 March 1982 (warnings 45 - 51, 
Table 17) during the looping phase. 

In contrast, Typhoon Marge is a case in which there are 


larger warning position errors for the objective technique 
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and smaller 24-h forecast errors from those warning posi- 
pons (Table 18). Objectively initiated  CLIPER forecasts 
performed better from 12 GMT 31 October through 12 GMT 02 
November 1983 (warnings 5 - 13, Table 13) while Marge was in 
the formative stages (< 65 kt) and then again from 13 GMT 03 
November through 00 GMT 05 November 1983 (warnings 18 - 23, 
Table 18) during the initial stages of recurvature. In both 
periods, Typhoon Marge was about to make a major course 
change. In this case, the JTWC warning positions fell on the 
wrong side of the turn as , shown previously in Fig. 3.7 and 
the corresponding CLIPER track departs significantly from 
the actual storm movement. 

The characteristics of Tropical Storm Winona (1982) and 
Typhoon Dot (1982) are very similar in both time of occur- 
rence (one month apart) and storm track (Typhoon Dot track 
approximately 5° north of Tropical Storm Winona). In both 
storms, the average initial position error evaluation 
resulted in a tie between the objective technique and the 
JTWC procedure (Winona, 22 n mi; Dot, 20 n mi). In the 24-h 
forecast evaluation, JTWC's warning positions resulted in 
smaller forecast error for Tropical Storm Winona ( 99 n mi 
compared to 110 n mi for the objective scheme). However, 
the objective technique's initial positions provided a 
superior forecast position error for Typhoon Dot ( 93 n mi 
compared to 105 n mi for forecasts from the JTWC warning 


positions). 


O SENSITIVITY TO STORM INTENSITY 


The objective warning position errors or the corre- 
sponding 24-h forecast errors are examined for Sensitivity 
to storm intensity. The 30-storm data base is stratified 
into storms of >130 kt (Super Typhoon), >65 kt (Typhoon) and 
<65 kt (Tropical Storm and Depressions). It should be noted 
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that the storms were of the indicated intensity for only 
part of their lifetimes, which may affect the results. 
Table 19 lists the total number’ of storms in each stratifi- 
cation and the percentage of win, tie or loss for both 
average initial position error and 24-h forecast position 
error. The most intense storms had the lowest percentage of 
win and tie cases for the objective scheme for both evalua- 
607 This is attributed to the most intense storms are 
also the easiest to locate by the JTWC and are given 
increased emphasis on their track forecasts. The objective 
technique may well be performing at the same level but would 
appear to be degraded in relation to JTWC in these intense 
Storm situations. Table 19 also indicates that good warning 
positions are produced by the objective technique for storms 


of Typhoon intensity and lower. 
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TABLE IO 
EE Z4-h forecast error category 
it 


(win, e, loss) versus initial position 
٭‎ ۰۶3 ۶۰۰۰۰۹ 66 ۳٣ 
24-h Forecast Category 
WIN DTIC TOSS 
initial WIN 7 2 6 
Position TIE 2 0 3 
Category LOSS 2 0 7 
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ա Example of a larger forecast error 
from a warning Ros Pon ا‎ ) that is more 
accurate than the official position (B). 
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TABLE LÎ 


Summary of warning position errors (n.mi.) for 
JTWC (JT) and objective (OBJ) for the independent 


sample storms during 1981. Entries are similar 
to Table 9. 
Storm No. of Avg. Std, Dev, wis T 
Warnings JT OBJ JT OBJ 
TY June 2 21 16 14 3 W -1Լ. 18 
Ts Roy 19 22 22 15 16 1 .06 
Ty Claira - 34 22 16 19 8 W -1.67 
STY Elsie D IES ES IET 1 W -2.00 
TY Hazen 40 24 20 13 و‎ ው ፤፡ ሽሻ... 
TS Jeff 20 34 29 20 2 W -. 88 
RIE 45 17 16 12 | ሃ -.38 
Total 214 | 
Weighted Average 22 18 14 11 
TABLE IZ 

Summary of warning position errors (n.mi. ք 

JTWC (JT) and objectives gen Lor the Praen E 

sample storms during 189821 ntries are similar 

to Table 9 
Storm No. of Avg Std. Dev. WTL T 
Warnings JT OBJ IT OBJ 

TY Nelson S 23 17 16 10 Wus += 2 que 
TS Winona 20 22 22 112 13 T Os 
STY Bess 45 16 2 10 1. Ն 1.43 
TY Judy 30 19 18 27 14 Ւ... 
TY Ken 36 14 15 6 2 Լ, 2 
TY Nanc E 12 11 11 6 W ON 
TY Pamela 68 ጋመ መጩ በከ Too o 
Total 285 
Weighted Average 19 18 15 کا‎ 


* Difference between the objective error and JTWC 
"٣*٢ error is significant at 95% confidence 
eve 


48 


م 


እን ct warning 


JTWC ( i and objective (OBJ 
sample storms during 1983. 
to Table 9 
Storm Noa On Avg. 
Warnings ST OBJ 
MY Tip 17 13 12 
TS Ben 112 17 34 
STY Forrest SL 10 15 
TS Georgia 13 10 19 
TY Lex Iu 18 18 
TY Percy 24 22 27 
TS Sperry 9 34 30 
TS Thelma 15 32 25 
Total 138 
Weighted Average 13 21 
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* Difference between the objective error and JTWC 
© cial error is significant at 93% confidence 


level 


JTWC (JT) 


positions for the 198 


Storm No. of 

Warnings 
TY June 6 
TS Ro 2 
mx Clara 25 
STY Elsie 25 
TY Hazen 33 
TS Jeff 10 
EY Kit SL 
Total 52 


Weighted Average 


TABLE 14 


Summary of 24-h forecast errors 
created with best track (BT) position, 


and objective 
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TABLE T5 


summary of 24-h forecast errors (m mi oom 
CLIPER initiated with best track (BT) positions, 
JTWC (JT) and objective (OBJ) warning 
positions tor thesrysezZ= scorns - 


Storm Nog GE Avg. Std. Dev. WTL T 
Warnings BT JT OBJ BT JT EORI 

TY Nelson 46 102 111 124 50 65 54 1. 1.02 
TS Winona 11 80 و‎ ۵ 52 DUDO E . 30 
STY Bess 38 1125 աա 61 65 68 ٢ ¿208 
T Judy 25 IA 96 63 64 69 W -.30 
TY Ken 28 71 87 86 49 و‎ ٢ W — 09 
TM ን 24 ER 66 66 20 355808 T -. QUE 
TY Pamela 56 122 134 143 աա:  .-: Ն ጋው 
"Eo 220 

Weighted Average 96 ԱԱ ԱՏ 58 64 69 

TABLE 16 


Summary of 24-h forecast errors (n.mi.) from 
CLIPER initiated with best ር ስ ን፡ ፓፕ bR NIPDOsItISmN 
JTWC (JT) and objective (OBJ) warning 
positions for the 1983 storms. 


Storm Noc ٤ Ave. Sid. Dev. WTL T 
c Warnings BT JT OBJ BT JT OBJ 
TY Tip 0 41 58 39 ጋን E Oo JE .07 
TS Ben 5 254 260 ጋያ) Wess INN 0 i , 1 
STY Forrest 24 74 79 89 46 54 54 E -09 
TS Georgia 6 63 58 89 20 13 46 E 56 
TY Lex 10 77 75 98 ከለ چ0‎ L 75 
TY Percy 17 192 2215. (55 50 47 53 ն es 
TS Thelma 6 135 141 143 104 112 87 E . 04 
Total 78 
Weighted Average lll 120 134 56 60 66 
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TABLE 17 
h forecast position errors 


statistics (n mi) for Typhoon Nelson 


y of 24 


Summar 


Obj-JTWC 


Wrng. 


Obj. 


Best TracK JTWC Wrng. 


Wrng. No. 
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TABLE 18 


Summary of 24-h forecast position errors 
statistics (n mi) for Typhoon Marge 


Wrng. No. Best Track JTWC Wrng. Obj. Wrng. Obj-JTWC 


E 49 128 93 35 
6 65 41 1 -30 
7 152 SH 88 -224 
8 172 205 215 10 
9 108 122 195 70 
10 39 187 32 E 
1 1 74 87 74 SES 
1.2 71 147 71 -76 
ÉS 1:59 198 143 -56 
14 183 258 238 0 
15 139 . 108 1023 50 
16 96 76 136 60 
17 61 31 67 36 
18 10 32 17 ze 
19 99 122 62 - 60 
20 114 147 137 -10 
21 91 77 106 2 
7 8l 134 112 > 
23 80 134 90 -44 
24 187 240 231 =, 
2.9 249 247 217 30 
26 360 د‎ 1 406 29 
27 546 487 545 58 
Avg. Disp. 1997 170 155 
Std. Dev. 117 112 12% 
TABLE 19 
Performance of objective technique 
based on intensity stratification 

Storm Intensity No. of Warning Forecast 

Storms Position Position 

> 130 Kts 5 20 80 20 80 

> 65 Kts wen 82 18 53 47 

< 65 Kts 8 75 25 43 57 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 


An objective technique for determining the warning posi- 
tion of a tropical cyclone has been developed and shown to 
be a viable forecasting "tool" for the tropical cyclone 
forecaster. The thought processes of the TDO have been 
Synthesized by the objective method of generating the 
"working best track" and future 24-h storm track, and the 
inclusion of the spatial and temporal weighting factors. 
This technique provides a consistent "first-guess'" for the 


inexperienced forecaster (Morford, 1979) or a first stepiin 


an objective format (Simpson, 1971; Elsberry, 1934) for 
forecasting tropical cyclone movement. Compared to the JTWC 
warnings, the three year (1981 - 1983) sample indicates 


slightly more accurate and consistent warning positions from 
the objective scheme. On a storm-by-storm basis, 21 of the 
30 storms had either improved or comparable warning position 
errors. 

JTWC's method of blending the +6-h forecast positions 
from the objective aids with the latest fix information to 
determine the warning position, which was initiated in 1933 
(Sandgathe, 1985), is similar in principle to this objective 
technique. During 1983, JTWC's method resulted in a smaller 
average warning position error. However, the objective 
technique still provided a more consistent position as indi- 
cated by the smaller standard deviation. 

Nine of the 30 storms tested had larger warning position 
errors than JTWC. Of these nine, five storms either looped 
or departed significantly from climatological tracks. It is 


these difficult forecasting scenarios in which the TDO 
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requires the best possible guidance. The CLIPER method does 
not include synoptic fields as predictors (Neumann, 1972) 
and thus will require additional guidance to handle these 
situations. 

The overall evaluations based on the 24-h forecast 
comparisons were not as promising. The total average fore- 
cast error was slightly smaller than JTWC's when using the 
objective scheme, due mainly to the reduced errors during 
19817 Only 45% of the 30 storms had improved or comparable 
forecast errors when CLIPER was initiated with the objective 
warning position instead of the JTWC warning position. Te 
is difficult to explain this diminished performance as exem- 
plified by the tracks of Tropical Storm Winona (1982) and 
Typhoon Dot (1982) described in section III-B. Comparing 
the standard deviations of the 24-h forecasts, it s 
apparent that there is not a significant difference in the 
consistency of the forecast positions in the 30 storms 
tested. Interestingly, the JTWC average 24-h forecast error 
would have been reduced in 16 of the 30 storms if the CLIPER 
forecast from the warning position had been used as 
guidance. 

The objective technique when applied to the strongest 
storms (> 130 kts) had the poorest performance relative to 
JTWC in terms of both warning positions and 24-forecast 
errors. This may be attributable to intense storms normally 
being characterized by well-defined eyes and circulation 
patterns that are easier to define accurately. JTWC will 
also place more emphasis on their forecasts for storms with 
more destructive capability. It may not be that the objec- 
tive technique is performing worse, and that JTWC's perform- 
ance is improved for intense storms. | 

The objective technique will require subjective enhance- 
ment by the TDO in looping scenarios and during periods of 


rapid acceleration after recurvature. Finally, "the 18ا0‎ 
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model is readily compatible to desk-top computers ոմ 
provides a fast, interactive forecasting model for the TDO. 
Optimally, the TDO will be able to quickly generate and 
evaluate a revised "working best track", warning position 
and short-term forecast track every time new fix information 


lS received. 


B. RECOMMENDATIONS 


The CLIPER regression coefficients were derived fora 
limited time and space domain. New coefficients should be 
derived for a sample that includes JTWC's entire area of 
responsibility. CLIPER contains no synoptic field informa- 
pon or physical interpretations. Dynamic models such as 
the One-way (Interactive) Tropical Cyclone Model (OTCM) or 
the Nested Tropical Cyclone Model (NTCM) incorporate both 
Synoptic data and physics and I would recommend that the 
objective scheme be coupled with each of these dynamic 
models to determine which pairing results in the greatest 
reduction of both warning position error and forecast posi- 
tion error. Finally, a real-time study of the objective 
scheme utilizing fix information and official JTWC warning 


positions should be conducted as a follow-on to this study. 
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